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Universal Multifractality in Quantum Hall Systems with Long-Range Disorder 

Potential * 



Rochus Klesse and Marcus Metzler 
Institut fur Theoretische Physik der Universitat zu Koln 
D-50937 Koln, Germany 

We investigate numerically the localization-delocalization transition in quantum Hall systems with 
long-range disorder potential with respect to multifractal properties. Wavefunctions at the tran- 
sition energy are obtained within the framework of the generalized Chalker-Coddington network 
model. We determine the critical exponent «o characterizing the scaling behavior of the local order 
parameter for systems with potential correlation length d up to 12 magnetic lengths /. Our results 
show that «o does not depend on the ratio djl. With increasing d/l, effects due to classical perco- 
lation only cause an increase of the microscopic length scale, whereas the critical behavior on larger 
scales remains unchanged. This proves that systems with long-range disorder belong to the same 
universality class as those with short-range disorder. 



PACS: 71.30.+h, 73.40. Hm, 71.50. +t 

Since the discovery of the Quantum Hall Effect (QHE), 
the localization-delocalization (LD) problem in two- 
dimensional (2D) disordered electron systems in a strong 
magnetic field has attracted much attention. It is now 
generally accepted that in the case of spinless non- 
interacting electrons only states in the center of the dis- 
order broadened Landau bands are extended, while at all 
other energies the states are localized [1]. The role of frac- 
tality in this transition was already observed by Aoki [2] 
investigating inverse participation numbers. The inher- 
ent connection between localization transitions and mul- 
tifractality was first recognized by Castelani and Peliti 
[3]. They have calculated generalized inverse participa- 
tion numbers (P(q) = J \ip(r)\ 2q d 2 r) and found a non 
trivial (/-dependence of the scaling behavior from which 
they concluded that wavefunctions at the critical point 
exhibit multifractal properties. Since then the multifrac- 
tal analysis of the LD transition has been developed con- 
siderably [4]. In the case of the QHE, Pook, Janssen [5] 
and Huckestein, Schweitzer [6] have calculated the /(a)- 
spectra of critical states and found them to be universal 
for quantum Hall systems with several types of short- 
range disorder potentials. Systems with long-range dis- 
order, where the potential correlation length d is large 
compared to the magnetic length I, have not yet been 
investigated, since the numerical effort increases rapidly 
with the ratio d/l within the framework of the models 
that have been used up to now. 

Indeed, the case of d I is especially interesting, be- 
cause then a very appealing semi-classical picture applies 
in which the LD transition is connected to the 2D per- 
colation transition [7]: In a smooth random potential V 
and strong magnetic field a classical 2D-electron executes 
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a fast cyclotron motion on a circle of radius I around a 
guiding center, which drifts along a contour V(r) = E'. 
For this reason the probability density of eigenstates with 
energy E = E' + hto c /2 is concentrated on strips of width 
I along these contours. An overlap between strips of 
separate contours at saddle points allows the electron 
to tunnel between them. Since this overlap decreases 
with decreasing I, tunneling becomes completely negligi- 
ble in the limit d/l —± oo. In this case the eigenstates 
are localized on single contours V(r) = E' with criti- 
cal size £ p diverging cx \E' — E c \ Vp as E' approaches the 
percolation threshold E c . Here v p = 4/3 is the criti- 
cal exponent belonging to the 2D-percolation transition. 
However, the LD-transition is - at least for small ratios 
d/l - characterized by a correlation length diverging with 
a different critical exponent v = 2.3 [1]. So the question 
arises, whether at large (but finite) ratios d/l effects due 
to classical percolation change the nature of the quantum 
mechanical LD-transition or not. 

We investigate this problem by analyzing critical states 
of generalized Chalker-Coddington networks with respect 
to multifractality. These models have been introduced 
by Chalker, Coddington [9] and generalized by Chalker, 
Eastmond [10] and Lee et al. [11] to describe quantum 
Hall systems with long-range random potential and are 
therefore suitable for our purposes. We analyze critical 
wavefunctions of network models corresponding to ratios 
d/l up to 12 and find that the characteristic multifractal 
exponent ctQ is independent of djl and the same as for sys- 
tems with short-range disorder investigated in [5,6]. So 
our numerical results directly confirm systems with d/l 
varying from up to 12 to be in the same universality 
class. Furthermore, a discussion of the influence of clas- 
sical percolation effects leads us to the conclusion that 
this remains valid even for arbitrarily large (but finite) 
ratio d/l, if the system size is sufficiently large. 

The network model consists of 2 x 2 unitary scattering 



1 



matrices representing saddle points of a smooth random 
potential which are arranged on a square lattice and con- 
nected by one-dimensional unidirectional channels, called 
links, corresponding to equipotential lines between the 
saddle points. The scattering matrices Sj describe the 
transitions from electron states of incoming links tpk^h 
to outgoing link states 4'm,ipn (see fig 1), 

(^ r ) = Sj (^), s j = ( Y \ ml ) . (i) 

\ Vn J \1PI J \t„k t ni J 

Traversing a link k an electron acquires an additional 
phase <pk, randomly distributed but constant for each 
link, which we absorb in the complex scattering coef- 
ficients t m k,t n k- The transmission amplitudes Tj = 
\t m k\ 2 = \t n l\ 2 , Rj = \t m l\ 2 = \t n k\ 2 are determined by 
the difference between electron energy E and random 
saddle point energy Uj G [— W, W], T^ 1 = 1 + exp((wj — 
E)/A),Rj = l — Tj. The relevant parameter is the ra- 
tio of the width W of the saddle point distribution and 
the tunneling energy A; W/A = yields the original 
Chalker-Coddington model, whereas W/A —± oo corre- 
sponds to classical percolation, since then the transmis- 
sion amplitude is zero or unity at every saddle point. 

To compare this model with explicit smooth disorder 
models we have to relate W and A to the magnetic length 
I and to the strength Vq = (V 2 ) 1 / 2 and correlation length 
d of the disorder potential. The energy distribution of 
the saddle points is limited by the potential strength Vq, 
therefore W « Vq. For a real saddle point of the ran- 
dom potential with curvature c = \V xx V yy \^l 2 the inverse 
transmission amplitude is 1 + exp((uj — E) 2it/cl 2 ) [12], 
so with the estimate c « Vo/d 2 one gets A « Vo(l/d) 2 
and W/A « (d/l) 2 . 

A wavefunction on the network is given by its com- 
plex amplitudes ipi on each link I, i. e. ^ G C N , where TV 
is the number of links. It is stationary at energy E if the 
scattering condition 

■0m = t m k1pk + tmirfu (2) 

is satisfied at each saddle point. For our purpose it is con- 
venient to use the disorder and energy dependent unitary 
operator U defined by 

Uei = t m ie m + t n ie n (3) 

(see fig. 1), where ej denotes the unity vector whose com- 
ponents vanish everywhere except at link j [13]. Now 
condition (2) simply states 

= (4) 

For a given disorder configuration this equation has so- 
lutions only for discrete energies. Instead of explicitly 
looking for these energies, which would be very cumber- 
some, we fix E and solve the more general eigenvalue 
equation Ue^ = e lv< &- A solution \I> v of this equation 



is obviously a stationary wavefunction at energy E for 
a disorder configuration corresponding to U' := e~ zlp U , 
i.e. a configuration that deviates from the initial one only 
by a simple global shift of the random phases. The lat- 
ter transformation does not influence the randomness of 
the coefficients, hence it does not change the statistical 
properties of the wavefunctions we are interested in. 

The square amplitudes of the normalized wavefunc- 
tions can be subjected to a multifractal analysis [4]: di- 
vide the system of size L into N = A -2 boxes bi(X) of size 
Lb = XL, X £ [0, 1] and determine the scaling behavior 
of the averaged q-moments of the box-probabilities 

( P9 ( X )) = ^J1 P ^ P *= I d 2 rMr)\ 2 , (5) 

with respect to A for real values of q. The scaling behavior 
can be described by a function r(q) that is defined by 

(P q (X)) <x X D+T ( q \ (6) 

where in our case D = 2. r(q) obviously increases 
monotonically with q due to normalization and has - less 
obviously - negative curvature, see e.g. [4]. Therefore its 
Legendre transformed, the so called /(a)-spectrum, 

dr 

a(q) = —, f(a(q)) = a(q)-q-T(q). (7) 

is a single-humped, positive valued function, restricted 
to D+oo < q < -D_oo, where -D±oo = &(<} ±oo). The 
/(a) -spectrum describes the scaling behavior of the en- 
tire distribution function of box probabilities and takes 
its maximum value D at ctQ = a(q = 0). This value ao 
is the critical exponent describing the typical scaling of 
the box-probabilities, Pt yp = exp (InP) c< A" . In phase- 
coherent systems the distribution of the box probabilities 
is of /og-normal type corresponding to a parabolic /(a)- 
spectrum [4] 

Indeed, this parabolic approximation, fitted only by ao, 
describes the /(a)-spectrum obtained in [5,6] in a broad 
range around the maximum. For these reasons ao is the 
most important scaling exponent and least influenced by 
statistical fluctuations. Other characteristic exponents 
are -D±oo describing the scaling of the minimum (-) and 
maximum (+) box-probabilities, P± cx A D± °°, and the 
exponent r(2) of the inverse participation number. The 
latter is connected to the anomalous diffusion exponent 
rj [14] via the scaling relation rj = D — r(2) [5,8]. 

We analyze wavefunctions at the critical energy E = 
for A = 1 and W varying from to 200 for networks of 
150 X 150 saddle points on a cylinder. As expected, the 
amplitude distributions of these wavefunctions are of log- 
normal type and a multifractal analysis can be applied. 
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For W < 75 we obtain almost the same /(a)-spectra 
as for systems with short-range disorder potential (see 
fig. 2). For W > 75 the sample to sample fluctuations and 
hence the statistical error increases rapidly, especially for 
large values of \q\. In order to obtain proper /(a)-spectra 
we had to choose larger minimal box sizes for increasing 
values of W. At W « 200 even for q = this minimum 
box size approaches the system size, and thus for W ^ 
150 no /(a)-spectra can be obtained. 

The exponent a shows no significant dependence on 
the width W, even for the largest values (see fig. 2). 
The value averaged over all systems is ctQ = 2.27± 0.01. 
For the original Chalker-Coddington model (W = 0) 
we obtain ao = 2.26 ± 0.01. This is in good agree- 
ment with the values 2.30 ± 0.07 and 2.29 ± 0.02 ob- 
tained for critical states in short ranged disorder sys- 
tems in [5] and [6], respectively. Furthermore, we find 
D_oo = 3.6 ± 0.2, = 1.0 ± 0.2 and the anomalous 
diffusion exponent r\ = 0.4 ±0.1 also very close to those 
in [5] (D-oo = 3.7 ± 0.1, = 0.95 ± 0.1) and [14,15] 
(r) = 0.38 ± 0.04). These results confirm the hypothesis 
of universal multifractality for long-range disorder poten- 
tials with correlation lengths up to d « -\/l50l « 121. In 
the following we argue that this scaling behaviour does 
not change in the limit L — t oo even for arbitrarily large 
ratios W/A and d/l. 

With this end in mind we discuss the influence of 
classical percolation on this LD-transition. The original 
Chalker-Coddington model (W = 0) exhibits no classi- 
cal percolation, because all saddle points are at the same 
energy, whereas in the generalized model percolation is 
incorporated and defines a new length scale. The micro- 
scopic length of the original model is simply given by the 
lattice constant a corresponding to the potential corre- 
lation length d. For the generalized model with W > A 
the microscopic length is significantly enhanced by the 
occurrence of saddle points whose energies Uj differ more 
than A from the electron energy E = 0. At these sad- 
dle points the coefficient Rj (or Tj) is exponentially small 
(~ exp(— |uj|/A)) and the amplitude in an incoming link 
is essentially transmitted in a single outgoing link only 
(classical behavior). Thus the amplitude of the wave- 
function remains constant on percolating paths avoiding 
saddle points with energy \uj \ < A. The typical distance 
covered by these classical paths sets a new microscopic 
length scale as a function of W/A, a(W/A). Figure 2 
clearly shows the appearance of the new length scale: the 
areas of approximately constant amplitude are consider- 
ably larger in the grey-scale plot on the right than on 
the left corresponding to wavefunctions of systems with 
W/A = 150 and W/A = 0, respectively. Identifying 
these classical paths with links, a model with arbitrarily 
large W can be rescaled to a model with lattice constant 
a' = a(W/A) and W' « A [16]. So the influence of clas- 
sical percolation can be eliminated by renormalizing the 
lattice constant and universal behavior is to be expected 



as long as the system size is larger than a(W / A). This 
explains the above-mentioned necessity to choose larger 
minimum box sizes, since they have to be larger than the 
microscopic length, and the breakdown of the multifrac- 
tal analysis when a(W / A) reaches the system size. 

To summarize our findings, we gave numerical evidence 
that systems with potential correlation lengths up to 12 
times the magnetic length have the same critical expo- 
nent ao as short-range disorder systems. We also argued 
that a will not change even for much larger ratios of 
d/l as long as the systems are sufficiently large. This is 
based on the observation that classical percolation which 
could be responsible for such a change only leads to an 
increase of the microscopical length scale that can be ac- 
counted for by rescaling the system. This leads us to the 
conclusion that the LD transition in quantum Hall sys- 
tems is universal for disorder potentials with arbitrary 
correlation length. 
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B. Huckestein, M. Janssen, D. G. Polyakov and M. Zirn- 
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Forschungsgemeinschaft for supporting this work within 
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FIG. 1. The Chalker-Coddington network. At each saddle 
point the scattering matrix Sj describes the transition from 
incoming to outgoing states. The operator U maps each in- 
coming link amplitude to the two outgoing links. 

FIG. 2. The ao values for different saddle point distribu- 
tions. (W is plotted logarithmically.) The insets show the 
/(a)-spectra of the wavefunctions and their amplitudes in a 
grey-scale plot, where darker areas denote lower square am- 
plitude. On the left one can see spectra for W = to 5 
compared to the parabolic approximation and one wavefunc- 
tion for W = and on the right W = 30 to 75 with W = 2 as 
reference and the plot depicts a wavefunction for W = 150. 



